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Abstract: The experimental technique of fluorescence light observation is used in current and planned 
air shower experiments that aim at understanding the origin of ultra-high energy cosmic rays. In the 
fluorescence technique, the geometry of the shower is reconstructed based on the correlation between 
viewing angle and arrival time of the signals detected by the telescope. The signals are compared to those 
expected for different shower geometries and the best-fit geometry is determined. The calculation of the 
expected signals is usually based on a relatively simple function which is motivated by basic geometrical 
considerations. This function is based on certain assumptions on the processes of light emission and 
propagation through the atmosphere. For instance, the fluorescence light is assumed to propagate with 
vacuum speed of light. We investigate the validity of these assumptions and provide corrections that can 
be used in the geometry reconstruction. The impact on reconstruction parameters is studied. The results 
are also relevant for hybrid observations where the shower is registered simultaneously by fluorescence 
and surface detectors. 



Introduction 

Since the 60's, when the use of extensive air 
shower (EAS) fluorescence light yield for the ul- 
tra high energy cosmic rays (UHECR's) detection 
was first proposed, many past, current and future 
experiments [1, 2, 3] utilize the effect to get a clue 
about the origin of cosmic rays. 
First applied for the Utah Fly's Eye detector [4] the 
emitted fluorescence light of EAS is used to recon- 
struct the shower geometry. The standard single- 
eye fitting procedure for the shower core location 
and direction starts with the determination of the 
plane containing the shower axis and the center 
of the eye (cf. Fig. 1). This so called shower- 
detector plane (SDP) fit uses tube pointing direc- 
tions, together with signal integrals. To determine 
the shower orientation within the SDP a correla- 
tion between viewing angle x% ar, d firing times of 
the tubes are used and a best-fit geometry with the 
expected arrival time at the telescope is accom- 
plished. The calculation of the expected light ar- 
rival time is motivated by basic geometrical con- 
siderations. 

Let t™ 1 be the time of light emission at point Si on 
the shower axis. Then, tf m is calculated as 




Figure 1 : Shower/Detector Geometry 
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where to is the time at which the shower passes 
the closest point at distance R p to the detector and 
Oi is the angle in the SDP which the shower axis 
makes with the i th pixel viewing towards Si (cf. 
Fig. 1). The shower front is assumed here to prop- 
agate with the speed of light in vacuum c. The time 
ti when the light reaches the telescope at pixel i 
viewing towards Si is then given by: 
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Where T f rop is the time of propagation from Si to- 
wards the telescope. Assuming also a propagation 
speed of cl- L = c for fluorescence light we obtain 
with (1) and (2) the "classical" formula used so far 
[5], 
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where xo and Xi are angles within the SDP, see 
Fig. 1. The best fit parameters R p , to and \o m 
Equation (3) are then found by minimizing a x 2 - 
function. The uncertainty of the three parameters 
depends on the particular shower geometry and is 
propagated also for the determination of primary 
energy and depth of shower maximum. In Eqn. 
(3), it is assumed that everything 

• propagates with vacuum speed of light, 

• takes place instantaneously 

• propagates on straight lines. 

We check the validity of these assumptions and 
provide corrections to Eqn. (3). 

Reduced speed of light 

Using high-speed electronics with trigger times in 
the scope of ns, the derivation of the expected light 
arrival time by basic geometrical considerations 
has to be revisited. The propagation speed of flu- 
orescence light v = c/n is reduced by an index 
of refraction n > 1 which is again a function of 
the traversed medium and wavelength A. The main 
emission lines of fluorescence light cover a wave- 
length range between 300 and 400 nm. Within that 
interval a change of n(A) is < 3% and negligible. 
Knowing the local density in air we can calculate 
the index of refraction n as follows [6]: 
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Figure 2: Arrival time difference \t rea i — t 
The telescope is located at position [0,0] at 1416 m 
a.s.l. (corresponding to the altitude of the Pierre 
Auger Observatory). The straight lines indicate the 
field of view between 1 and 3 1 deg 
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where h is the height a.s.l and bj as well as Cj at- 
mospheric parameters for the U.S. standard atmo- 
sphere given in [6]. 

To estimate the impact of a realistic atmosphere we 
calculated the difference of the light arrival times 
between the cases of vacuum and realistic speed 
of light from different parts of the atmosphere as 
shown in Fig. 2. We can see that the time differ- 
ence is larger for light propagating near the earth 
surface as we would expect and that differences of 
more than 20-25 ns can occur. Furthermore, we see 
that coming in showers, with respect to the tele- 
scopes, are expected to have a more constant offset 
contrary to going away showers. 

De-excitation lifetimes 

Another additional delay to the expected arrival 
time are excitation and de-excitation processes 
within the shower development induced by low 
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Figure 3: Lifetime of individual transitions as a 
function of height a.s.l. for dry air. The thick 
line indicates the weighted lifetime according to 
different intensity fractions. The width of that 
line roughly corresponds to a ±40 K temperature 
change. 



With Eqn. (1), (2) and (6) the new expected arrival 
time ti can then be written in the form 



c- \sm(xo-Xi)J 
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d i denotes an averaged reduced speed of light for 
the particular path i, Eqn. (7) replaces the classical 
formula Eqn. (3) for fitting the shower geometry. 
Studies on simulated CORSIKA events taking into 
the Pierre Auger detector MC show, that differ- 
ences in the arrival time of 20-30 ns and more can 
occur. The fit parameters R p , to and %o are more 
affected for coming in shower and differences of 5 
m, 30 ns and 0.05 deg emerge, respectively. 



energy electrons and positrons (^40 MeV). Al- 
most all of the air fluorescence in the wavelength 
range between 300 nm and 400 nm originates from 
transitions of molecular nitrogen N 2 or molec- 
ular nitrogen ions which make ~78% of 
the air composition. Excitation times are around 
10 -15 ns and, thus, negligible for our purposes. 
De-excitation lifetimes can be of the order of 30- 
40 ns but are affected by quenching [7] . Here exci- 
tated states transfer their energy into rotations, vi- 
brations or translations of other molecules without 
emitting optical photons. As a consequence one 
has to introduce an additional radiationless deacti- 
vation term r q . The inverse mean lifetime r„ be- 
comes 



T u (p,T) To, T c „(p,T)' 

where 1 /tq„ is the sum over all constant transition 
probabilities. The expected de-excitation lifetimes 
of the main transitions as a function of altitude are 
shown in Fig. 3. 

The functional form of the weighted line in Fig. 3 
can be parameterized in good approximation by 

Tv ( h ) = -h/H T -i (6) 

a ■ e n / M + 1 

with tq v — 37.5 ns, H — 8005 m and a = 95. 



Bending of light 

In addition to the time delay, the light path also 
changes according to Fermats principle resulting 
in an abberation of the viewing angle Xi- Snell's 
law states that the ratio of the sines of the angles of 
incidence a rea i and refraction a b s is equal to the 
inverse ratio of the indices of refraction 

rii sin a rea i = n r sin a obs . (8) 

Consequently the telescope "detects" the light 
higher in the atmosphere than it was actually pro- 
duced. The angle difference is shown in Fig. 4. 

An angular difference of a a b s — a rea i ~ 0.05 deg 
implies a shifted observed emission point of about 
h = 30 m higher at a distance of d = 30 km. 
That could cause a delay of ~100 ns for the ex- 
pected impact time on ground. This is particularly 
important for so called hybrid observations where 
the same EAS is detected by the ground array and 
a fluorescence detector realized at the Pierre Auger 
Observatory. The resulting relative timing offset 
has to be taken into account for an accurate re- 
construction. In order to estimate the impact of 
bended and delayed fluorescence light on the time 
offset a toy model is used. This toy model simu- 
lates air showers in vacuum and a realistic atmo- 
sphere and determines the expected impact time 
on ground via the detected fluorescence light at the 
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Figure 4: (Top) Illustration of the path deviation 
and the resulting impact on the observed emission 
point. (Bottom) Arrival angle difference between 
direct and curved path for a realistic atmosphere. 
The different colors indicate the difference a b s — 
a real in deg. 

telescope. Other effects like light attenuation, tele- 
scope characteristics etc. are neglected. The ex- 
pected impact time for a realistic atmosphere mi- 
nus the impact time for a vacuum atmosphere (time 
offset) as a function of 9 is shown in Fig. 5. There 
is an increase for inclined showers mainly caused 
by bended fluorescence light. The rise at low 9 is 
caused by the absence of quenching effects at high 
altitudes. 

Conclusion 

We studied the impact and influences of realis- 
tic light propagation within shower reconstruction 
from fluorescence light measurements. Differ- 
ences in the arrival time of 20 - 30 ns and more can 
occur, which affects the geometry and, to a minor 
extend, the profile reconstruction. For certain ge- 
ometries the time offset is expected to exceed 100 
ns. A correction for the "classical" fit function is 
provided that can readily be applied to reconstruct 
fluorescence detector data. 
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Figure 5: Expected time offset as a function of 
9 for the toy model. The impact point is always 
26 km away with varying starting points along the 
maximum viewing angle, here 30°. 
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